Formation of regular structures in the process of phase separation 
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Phase separation under directional quenching has been studied in a Calm-Hilliard model. In 
distinct contrast to the disordered patterns which develop under a homogeneous quench periodic 
stripe patterns are generated behind the quench front. Their wavelength is uniquely defined by the 
velocity of the quench interface in a wide range. Numerical simulations match perfectly analytical 
results obtained in the limit of small and large velocities of the quench interface. Additional pe- 
riodic modulation of the quench interface may lead to cellular patterns. The quenching protocols 
analyzed in this paper are expected to be an effective tool in technological applications to design 
nanostructured materials. 

PACS numbers: 47.54.-r, 64.75.-g, 05.45.-a 



The dynamics of phase separation in multiphase sys- 
tems has been investigated intensively during the last 
decades (J, Q • A main paradigm are the spontaneously 
arising spatial concentration variations of a characteristic 
average domain size when a homogeneous binary mixture 
is quenched into the thermodynamically unstable region. 
With the progress of time coarsening takes place, i.e., 
the average domain size increases continuously. The or- 
der parameter, usually the relative concentration of the 
two species, obeys a conservation law in contrast to the 
large number of hydrodynamic pattern forming systems 
[3|| . For technological applications it is highly desirable 
to run a phase separation process in a controlled man- 
ner to create regular structures. This is important to 
design nanostructured materials and nanodevices in di- 
verse fields, ranging from bioactive patterns [J] to poly- 
mer electronics [5j . The arrangement and the size of the 
domains which form the functional elements have a cru- 
cial impact on the device performance, e.g., in photo- 
voltaics, LEDs, and electronic circuits made of polymer 
blends [5|, la, 0|- Previous attempts to manipulate the 
phase separation in binary mixtures by various external 
fields did not lead to a satisfactory control of the pat- 
tern morphology [3, H], [lfj, [II], [l2[ . A first breakthrough 
has been established recently by introducing a persistent 
spatially periodic temperature modulation in a model of 
the phase separation in binary mixtures [1 31 ] . In this case 
stripe patterns with the periodicity slaved to the exter- 
nally imposed one can be stabilized against coarsening 
above some critical modulation amplitude. 

In this paper a new effective mechanism to create pe- 
riodic stripe patterns by directional quenching will be 
presented. Their wavelength is uniquely selected by the 
velocity of quench interface. If in addition a spatially pe- 
riodic modulation of the quench interface is introduced 
also cellular patterns can be generated. 

The appropriate mean field description of phase sep- 
aration is commonly based on the generic Cahn-Hilliard 
(CH) model [3 (model B [l5|). It has been widely used 
to study the dynamics of phase separation processes in a 
large variety of systems, such as binary alloys, fluid mix- 
tures, and polymer blends [l|, Q- In the one-dimensional 



case the CH model is described by the following PDE 

d t u = d xx (-eu + u 3 -d xx u) , (1) 

where u(x,t) is a real order parameter, e.g., in a binary 
mixture the difference of concentration of one species 
from that at the critical point and e is the control pa- 
rameter. According to Eq. (pj the spatial average, (u), 
of u(x, t) is conserved. 

To keep the analysis most transparent it is sufficient 
to concentrate on the case of (u) — (the so-called 
critical quench). The homogeneous solution u — be- 
comes unstable for e > against linear perturbations 
~ e at + l ix w ith wavenumber q £ (0, y/e) and growth rate 
(7 = q 2 (e — q 2 ). The most unstable (fastest growing) 



mode is characterized by q„ 



'e/2 with a r , 



74. 



For e > a one-parameter family of stationary periodic 
solutions u p (x,k) of Eq. (Q]), the so-called soliton-latticc 
solutions, exists which can be expressed in terms of the 
Jacobian elliptic function sn as: 



u p {x,k) = — - sn(-,fc) , £ = 



1 + fc 2 



(2) 



The modulus k £ (0, 1) is related to the wavenumber q: 

t/2 
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2K(fc)£ 



K(t) = /^ 
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k 2 ^ 2 



(3) 



sm <p 



where K is the complete Jacobian elliptic integral of the 
first kind. It is known, however, that any periodic solu- 
tion of Eq. (p} is unstable against period doubling, i.e., 
against coarsening [161 ]. In the limit q <C q m the growth 
rate of the corresponding destabilizing mode is given as 
MM: 



e 2 l6eyLp(-2irq m /q)/(irq m /q) 



(4) 



Note that e can be scaled out in Eq. |T]) which is reflected 
in the e-dependence of u p , q, and a p . Coarsening becomes 
extremely slow and the solutions of type ([2]) persist for a 
long time At p rs 1/c p at small q (k — ► 1). This situation 



is favorable to the general goal of controlling phase sepa- 
ration to achieve long lived periodic structures. However, 
with just a homogeneous quench from negative to posi- 
tive values e in Eq. ([1]) periodic structures with a single q 
can never be obtained: Any random perturbations about 
u = will develop initially into a superposition of the 
fastest growing Fourier modes with average wavenumber 
(q) rs q m . In the next stage coarsening sets in; (q) de- 
creases continuously in time and follows a scaling law 

(<Z>~l/log*pIEi 

Consequently, the question arises, how to control the 
quenching process such that a quasi-stationary periodic 
solution of type ([2]) can be generated. It is well-known 
that temporal and/or spatial modulations of the control 
parameter serve as a powerful tool to influence the pat- 
tern selection processes (see, e.g., [18|, LL9|). One of the 
simplest cases is directional quenching where a jump of 
the control parameter is introduced and the arising in- 
terface is dragged with a constant velocity. This has mo- 
tivated us to analyze systematically the impact of the 
directional quenching in the generic CH model. Thus 
we have achieved a clear understanding of controlling 
the generation of periodic structures. In fact, directional 
quenching has been used in earlier numerical simulations 
of certain phase separation model and some kind of regu- 
lar patterns have been observed |20l . [2l| , to which we will 
return below. In the CH model (JTJ) directional quenching 
is realized by changing e from a negative value at x < x q 
to positive one for x > x q , i.e., dividing the system into 
a stable and unstable region. The quench interface (re- 
ferred to as QI in the sequel) at x q is moving in the 
laboratory frame with a velocity v, i.e., 



e(x,t) = 



—e 

+e 



x < —vt , 
x > —vt . 



(5) 



Numerical simulations of the Id CH model JT]) with 
the directional quenching ([5]) demonstrate that a peri- 
odic solution develops behind the QI in the unstable re- 
gion. Typical examples for large and small velocities v 
of the QI are shown in Fig. [1] For v above some crit- 
ical value, v* , the periodic solutions detaches from the 
moving QI and the wavelength of the solution becomes 
independent of v [Fig. HJa)]. In contrast for v < v* the 
solution remains attached to the QI with the wavelength 
uniquely determined by v. Decreasing v the solution de- 
velops into a periodic kink lattice (sharp changes between 
u = is/e) where new kinks are continuously generated 
at X = x q (t) = —vt [Fig. QJb)]. The period of the so- 
lution, 2A, turns out to be uniquely defined by veloc- 
ity of the QI and is shown in Fig. [51 For v — > one 
has A ~ 1/v whereas for v > v* one finds A = ir/q*. 
Although the periodic solutions far away from the mov- 
ing QI are in principle unstable against period doubling, 
the coarsening is extremely slow for patterns generated 
with q <C q m [see Eq. Q]. Thus the extension L p of the 
(quasi-ideal) periodic solution behind the QI can be esti- 
mated as L p = vAt p rs v/(Tp where cr p is the growth rate 
of the unstable period doubling mode given in Eq. ^ . 



(a) 



(b) 

FIG. 1: Solutions of Eq. © with the QI © for e = 1 in the 
comoving frame with the QI at x = 0: v — 2 > v* (a) and 
v — 0.02 <C v* (b). Only a part of the system of total length 
l x — 4096 is shown. 




FIG. 2: Kink separation length A = iv/q multiplied with the 
velocity v of the QI as a function of v for e = 1 (full circles with 
a solid line as a guide to the eye); v* = 1.622 from Eq. (0. 
The doted and the dashed lines correspond to Eq. and 
Eq. (|12p , respectively. 



The two limiting cases of large and small velocity v of 
the QI can be captured analytically. For large v we con- 
sider for instance the initial condition u = everywhere 
except a hump u > localized near x = 0. Then the time 
evolution of this initial perturbation is governed by the 
motion of wave fronts to the left and to the right with a 
well-defined velocity v* and wavenumber q*. These quan- 
tities can be calculated by a linear stability analysis of 
the leading edge of the front in the comoving frame. One 
arrives thus at the so-called marginal stability criteria 
GHH: 



dv{q ,v*) * Im[a(q ,v*)] 
= 0, Rc[a(q ,v )] = 0, q = — ,(6) 

where a(q,v) — q 2 (e — q 2 ) + iqv and go is a complex 
number (saddle point). Eqs. ([6]) lead to the solution 

^7 + o/o \Va 



y- 



i) 



,3/2 



= 3(V7 + 3)»/ 2 cl/a 

8^(^+2) 



(7) 



The phase velocity and the wavenumber of the propagat- 
ing periodic solutions obtained from the numerical simu- 
lations of Eq. ((T|) for v > v* , which do not depend on v, 
agree perfectly with v* and q* given by Eq. |J7J (Fig. [2]). 
In the opposite limit v — ► our starting point is a 
particular stationary solution of Eq. ([T]) for v = inter- 
polating between u = at x < and u — \fk at x > 



which is characterized by a sharp front at x ~ 0. If the 
QI according to Eq. © starts to move, the sharp front 
would initially follow. But since the spatial average (u) is 
conserved regions with u < have to be generated in the 
region x > x q which leads to the formation of kink lattice 
[Fig. QJb)]. Its formation can be understood in terms of 
a fast switching stage and slow pulling stage: first a new 
kink is generated in a short time at x w x q . During the 
slow stage this kink is pulled by the QI whereby its am- 
plitude and the distance to the next kink behind, Ao(t), 
increase until Ao exceeds some limiting value Xo t max and 
then a new kink is generated (Fig. [3]). Repeating this 



ing expression 



Outer solution 




¥ 



(10) 



Starting from u — uq the outer solution u OM t will grow 
until at x = Xq (second kink) the maximal possible 
amplitude u max — y/e is reached. The maximal inter- 
val Xo t max — Xmax — Xmin where the outer solution is 
supported corresponds obviously to the minimal value 
Umin — \/2e/3 of «o [C = in Eq. (fT0|) ]. Let us now cal- 
culate the equilibrium distance between kinks A (e.g., the 
distance between the third and the fourth kink in Fig. [3]). 
During the slow pulling stage the distance between the 
second and the third kink increases until it reaches its 
maximal value A when a new kink is generated at the QI. 
Inspection of Fig. [3J shows that then because of (u) = 
the area Xy/e under the curve between two neighboring 
kinks away from the QI equals twice the area S max under 
the curve between the first and the second kink with 



FIG. 3: Sketch of the solution of Eq. JTJ with the QI © in 
the comoving frame for x > x q during the slow stage at the 
time shortly before the new kink will be generated at x q — 0. 
The domain of u ut(x) connecting to m n at x — xo is marked 
by the vertical dashes. 
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udx = 



s)~*- 



(11) 



With the use of the outer solution © one obtains from 

Ay£ = Aomax 



process a regular kink lattice develops in the wake of the 
QI with a kink separation length A (i.e., with the period 
2A), which is uniquely determined by the velocity v of 
the moving QI (Fig.©. 

During the slow stage the solution of Eq. (JTJ) in the 
interval < x < Ao(i) can be described in a comoving 
frame x — ► x + vt by 



vu 



d x {- 



eu - 



d x 



(8) 



which is obtained from Eq. (UJ by an x integration 
with the boundary conditions u — and zero flux 
d x (—eu + u 3 — d xx u) = for x = x q = 0. In addition 
the explicit time dependence dtu can be safely neglected 
since it is only relevant during the fast switching stage. 
The solution of Eq. ([5]) can be separated into a strongly 
varying inner solution Ui n (x) (0 < x < Xa, first kink) 
and to the almost flat outer solution u ou t(x) (x > xo) 
(Fig. [3j . For calculating m n the term vu sa in Eq. © 
and for u ou t the term d xx u sw can be neglected which 
leads to 



{d x Ui 



uj„ (e 



in x 

2 ' 



C + 2u in D 



C = {d x u in ) 2 \ x=0 , D = [d xx u m + u in {t 

Ox'Uout — o 9 

3<«t - e 



\X=XQ 7 

(9) 



The inner and outer solutions are matched at the point 
x = xq chosen such that "flatness" conditions d x ui n = 
dxxUm — hold. This gives for uq = Ui n (xo) the follow- 



9 v v 



(12) 



in perfect agreement with the results of numerical simu- 
lations in the limit v — > (Fig.©. 

The generalization of the analysis to the off-critical 
quench, (u) ^ 0, is straightforward. The expressions ([7]) 
for v* and q* hold except the replacement e — v e — 3(u) 2 . 
In the limit v — ► the distance A+ between two kinks 
in the region u > 0, and A_ for u < 0, respectively, 
become different. We find A + — A_ = (u)A/\/e and for 
the resulting period A of the kink lattice 



A = A+ + A_ = 



^ 



25f)(uf 



(13) 



As in the case (u) — 0, we found Eq. (|T3j) confirmed in 
numerical simulations of Eq. |T]) for different values of 
(u) in the limit v — > 0. 

Let us now switch to the two-dimensional case u(x, y, t) 
where we study first numerically the 2d version of CH 
equation |T]): 



d t u = V 2 (- 



(14) 



with the moving QI (|5|). Zero flux boundary conditions 
have been used at x = 0,l x and periodic boundary con- 
ditions at y = 0, l y . Initially the QI is located at x q = l x 
moving from right to left. The system size was l x = 512, 
ly = 256 and we start with the homogeneous solution 
u = (u) with small superimposed noise of the strength 
8u where 6u <C y/e and Su <C (u). Thus the well-known 



Ginzburg criterion, necessary for the validity of a mean- 
field description of a phase separation process [23| , is sat- 
isfied: in fact the dynamics does not depend on the par- 
ticular choice of Su. 

For the off-critical quench (u) ^ when v < v* always 
regular stripe patterns with domains parallel to the QI 
were found [see, e.g., Fig. Ufa)]. This situation is covered 
by an Id analysis presented before where the period of 
the structure is uniquely determined by the velocity of 
the QI. In the limit v — > the period of the patterns 
found in our numerical simulations agree with (|13[) . For 
v > v* irregular coarsening patterns similar to the case 
of a spatially homogeneous quench have been observed. 

In contrast, in the case of the critical quench (u) = 0, 
the orientation of the domains depends on the velocity of 
the QI. At small v periodic patterns with domains per- 
pendicular to the QI are formed [similar to Fig. lUJb)]. 
Then for v above v c w 0.45 the Id stripe patterns paral- 
lel to the QI appear as in the off-critical case described 
above. Finally v > v* leads eventually to irregular 
patterns similar to the case of a spatially homogeneous 
quench. 

The earlier studies of a related model subjected to di- 
rectional quenching in 2d [2CJ, [21| , which demonstrate the 
existence of regular patterns as well, were purely numeri- 
cal and do not give real insight into the pattern selection 
mechanisms. In particular, for the relevant off-critical 
case a noise strength Su > (m), violating the Ginzburg 
criterion, was chosen |2l| . Consequently, for instance, 
the change of the domain orientation from perpendicular 
to parallel with increasing v reported in [2l| is considered 
as an artefact. 

Finally we have studied the influence of a periodic 
modulation of the QI which reads as follows 



e(x,y,t) 



-e , x < l x + a cos(py) — vt 
-e , x > l x + a cos(py) — vt 



(15) 



In the case of a critical quench we found that the velocity 
v c at which the transition from parallel stripe patterns 
[similar to Fig. HJa)] to perpendicular ones [Fig. HJb)] 
occurs, depends on the modulation amplitude a. This 
dependence is very strong for values a of the order of the 
typical domain size at the initial stage of phase separa- 
tion (a ~ A TO = n/q m ). Furthermore, v c decreases with 
decreasing modulation wavenumber p. For p smaller than 
the wavenumber q m of the fastest growing mode patterns 
with a cellular morphology forming behind the moving QI 
have been found [Fig. IUc)]. In the case of an off-critical 
quench we found that (u) ^ favors the formation of 



regular cellular planforms [Fig. [4jd)] at intermediate QI 
velocities, in analogy to the transition from parallel to 
perpendicular stripes for (u) = 0. 

As a main result of the paper we have demonstrated 
that directional quenching in CH model leads to the 
formation of periodic solutions with the wavelength 
uniquely selected by the velocity of quench interface. 
This is in contrast to Ginzburg-Landau models with a 
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(b) 



(d) 



FIG. 4: Snapshots of the phase separation in 2d at the time 
when the QI almost reaches the left boundary. Straight QI 
© for e = 1, (u) = 0.1, v = 0.05 (a). Modulated QI (To) for 
e = 1, a = 4, p = tt/16: (u) = and v = 0.01 (b), v - 0.055 
(c); (u) =0.1, W = 0.1 (d). 



nonconserved order parameter where a moving quench 
interface leaves no kinks behind in the wake [24j . 

The Id CH model has recently found a new interest- 
ing application in liquid crystals to describe Ising walls 
between symmetry degenerated director configurations 
[2a . [26| . It should be possible to verify our predictions 
on the wavelength selection by dragging the liquid crystal 
layer into a homogeneous magnetic field. 

In conclusion, controlling phase separation by direc- 
tional quenching turns out to be a promising tool to cre- 
ate regular structures in material science. Although slow 
coarsening cannot be avoided by directional quenching in 
principle, long lived periodic patterns can be "frozen in" , 
e.g., by a deep quench, induced polymerization, chemi- 
cal treatment etc. In this respect it would be certainly 
rewarding to study in more depth the 2d case where in- 
teresting cellular structures have been detected and in 
particular the 3d case. 

I am deeply indebted to the late Lorenz Kramer from 
whom I have benefited a lot regarding the concept used 
in this paper. It is a pleasure to thank W. Pesch 
and W. Zimmermann for stimulating discussions and 
for critically reading the manuscript. Financial support 
by Deutsche Forschungsgemeinschaft Grant SFB 481 is 
gratefully acknowledged. 



[1] J. D. Gunton, M. San Miguel, and P. S. Sahni, in Phase 
Transitions and Critical Phenomena, edited by C. Domb 
and J. L. Lebowitz (Academic Press, London, 1983), Vol. 
8. 

[2] A. J. Bray, Adv. Phys. 43, 357 (1994). 



[3] M. C. Cross and P. C. Hohenberg, Rev. Mod. Phys. 65, 

851 (1993). 
[4] J. Voros, T. Blattler, and M. Textor, MRS Bull. 30, 202 

(2005). 
[5] H. Sirringhaus, Adv. Mater. 17, 2411 (2005). 



[6] D. C. Coffey and D. S. Ginger, J. Am. Chem. Soc. 127, [16 

4564 (2005). [17 

[7] G. Fichet, N. Corcoran, P. K. H. Ho, A. C. Arias, 

J. D. MacKenzie, W. T. S. Huck, and R. H. Friend, Adv. [18 

Mater. 16, 1908 (2004). [19 

[8] H. Tanaka and T. Sigehuzi, Phys. Rev. Lett. 75, 874 

(1995). [20 

[9] C. L. Emmott and A. J. Bray, Phys. Rev. E 54, 4568 [21 

(1996). 
[10] L. Berthier, Phys. Rev. E 63, 051503 (2001). [22 

[11] A. A. Golovin, A. A. Ncpomnyashchy, S. H. Davis, and 

M. A. Zaks, Phys. Rev. Lett. 86, 1550 (2001). [23 

[12] A. Voit, A. Krekhov, W. Enge, L. Kramer, and [24 

W. Kohler, Phys. Rev. Lett. 94, 214501 (2005). 
[13] A. P. Krekhov and L. Kramer, Phys. Rev. E 70, 061801 [25 

(2004). 
[14] J. W. Cahn and J. E. Hilliard, J. Chem. Phys. 28, 258 [26 

(1958); 31, 688 (1959). 
[15] P. C. Hohenberg and B. I. Halperin, Rev. Mod. Phys. 49, 

435 (1977). 



J. S. Langer, Ann. Phys. 65, 53 (1971). 

T. Kawakatsu and T. Munakata, Prog. Theor. Phys. 74, 

11 (1985). 

W. van Saarloos, Phys. Rep. 386, 29 (2003). 

S. Riidiger, E. M. Nicola, J. Casademunt, and L. Kramer, 

Phys. Rep. 447, 73 (2007). 

H. Furukawa, Physica A 180, 128 (1992). 

B. Liu, H. Zhang, and Y. Yang, J. Chem. Phys. 113, 719 
(2000). 

E. Ben-Jacob, H. Brand, G. Dee, L. Kramer, and 

J. S. Langer, Physica D 14, 348 (1985). 

K. Binder, J. Chem. Phys. 79, 6387 (1983). 

N. B. Kopnin and E. V. Thuneberg, Phys. Rev. Lett. 83, 

116 (1999). 

C. Chevallard, M. Clerc, P. Coullet, and J.- M. Gilli, 
Europhys. Lett. 58, 686 (2002). 

T. Nagaya and J. M. Gilli, Phys. Rev. Lett. 92, 145504 
(2004). 



